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Abstract. The three components of the heat-flux vector F = pC p {u'T') are numerically computed for a stratified rotating 
turbulent convection using the NIRVANA code in a flat box. The latitudinal component F s proves to be negative (positive) 
in the northern (southern) hemisphere so that the heat always flows towards the poles. As a surprise, the radial heat-flux F T 
peaks at the equator rather than at the poles (Taylor numbers O(10 6 )). The same behavior is observed for the radial turbulence 
intensity (u' r 2 ) which for free turbulence is also believed to peak at the poles (see Eq. 1 1 9i below). As we can show, however, 
the consequences of this unexpected result (also obtained by Kapyla, Korpi & Tuominen 2004) for the theory of differential 
rotation are small as mainly the F a is responsible to solve the 'Taylor number puzzle'. 

In all our simulations the azimuthal component F^ proves to be negative so that the rotating turbulence produces an westwards 
directed azimuthal heat-flux which should be observable. Fluctuations with higher temperature are expected to be anticorrelated 
with their own angular velocity fluctuations. We find this rotation-induced result as understandable as the F$ is closely related 
to the radial A-effect which is known to be also negative in stratified and rapidly rotating convection zones. 
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1. Motivation 

There are more and more data describing the surface rotation 
law of stars, the majority of which complies both with the sign 
and amplitude of the known equatorial acceleration of the Sun. 
In order to compare the various observational results a relation 



pole 



has been introduced with n" representing the ^-dependence 
of the equator-pole difference of the surface angular velocity. It 
proved to be weak in the first papers in this research (n" =t 0. 15, 
Hall 1991) while in recent studies values of 0.58 (Messina & 
Guinan 2003) and 0.66 (Reiners & Schmitt 2003) have been 
found. A rotational influence upon the rotation law is clearly 
existing but it seems not to be too strong. 

To find the influence of the global rotation on the turbu- 
lent convection in the stellar convection zones is thus the basic 
problem of the theory of differential rotation. The resulting flow 
pattern ('rotating turbulent convection') simultaneously trans- 
ports both the heat and the angular momentum. The resulting 
angular momentum transport has been described with the A- 
effect whose properties easily can be reproduced with box sim- 
ulations (see the references in Rudiger & Hollerbach 2004). In 
the present paper it is the eddy heat-flux in rotating turbulent 
convection which shall be rediscussed here by means of new 
box simulations. 

For the usual Boussinesq relation 



F = - P T X VS 



(2) 



between the eddy heat-flux and the entropy gradient Weiss 
(1965) and Durney & Roxburgh (1971) parameterized the in- 
fluence of a global rotation by 



(1) X =Xo 



+ ecos 0) 



(3) 
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with s representing the rotational influence. The eddy conduc- 
tivity is assumed to involve the main influence of the turbulent 
flow pattern which due to the basic rotation strongly depends 
on the colatitude 0. The ansatz leads to a ^-dependence of 
all thermodynamic quantities with a temperature difference of 
pole and equator so that a meridional flow is the immediate 
consequence. There was no possibility at that time for a de- 
tailed theory of the connection of xo and s with the character- 
istic parameters of the turbulence (see also Belvedere, Paterno 
& Stix 1980). 

In the magnetohydrodynamics of the Sun meridional flows 
are playing a more and more important role. The theory of the 
differential rotation in the convection zone and the tachocline 
as well as the theory of the solar dynamo are important ex- 
amples. We shall discuss details of the thermodynamics of the 
rotation theory in this paper. We shall derive the ideas for the 
structure of the eddy-conductivity heat tensor (next Section). 
In the following Sections the results will be compared with the 
results of nonlinear simulations with the finite-difference code 
NIRVANA. 
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2. Quasilinear theory 

In rotating turbulent fluids the relation between the turbulent 
heat-flux F = pC p (u'T') and the superadiabatic temperature 
gradient 

P=-L-VT (4) 
(g gravity) is a tensorial one, i.e. 

Ft = pC v XijPj (5) 

with p = Vr ad - VT = -(T/C p )V5 where S is the gas entropy. 
In the simplest case known from the mixing-length theory it is 
Xij - XT^ij so that Eq. results (instead of F = -pC P XT^T 
for incompressible fluids). 



to the poloidal components (0 of the heat-conductivity tensor. 
Hence, now 

Xrr = XT + (X + X\\)® 2 COS 2 8, 

Xer = ~X\\Q 2 sin 6 cos 8, (10) 

where some of the notations have slightly been changed. The 
sign of the unknown quantity x may fix the latitudinal profile 
of the radial heat-flux. 

In the following we shall compute JlQi for various rotation 
rates by means of numerical box simulations. We shall find, in- 
deed, that the pole-equator difference of the radial heat-flux is 
a sensitive function of rotation, stratification and/or boundary 
conditions. This is not true, however, for the latitudinal heat- 
flux which will prove to be towards the poles in all our simula- 
tions. 



2.1. Isotropic turbulence subject to rotation 

For rotating fluids without any other preferred direction than 
the rotation vector £2 it is simply 



Xij =Xr6ij +X\\ Q i Q ) +X £ ipi Q P 



(6) 



(Kitchatinov, Pipin & Riidiger 1994). The poloidal components 
of ijfji in spherical coordinates are 



Xrr =XT+COS 2 6x\\Q 2 , 

Xer = - sin #cos 8x\\Q 2 . 



(7) 



If the averages are taken over the horizontal plane, /3 has only a 
radial positive component which is positive in the convectively 
unstable zones. 

With the heat-transport in rotating but otherwise 
isotropic turbulence is rather simple. If the latitudinal heat- 
transport goes to the poles (equator) then the radial heat- 
transport at the poles is always stronger (weaker) than at the 
equator. In the first case the x\\ is positive and in the second 
case it is negative. 

The coefficients xt an d X\\ nave been computed by 
Kitchatinov, Pipin & Riidiger (1994, their Fig. 1). They both 
are positive. The^j vanishes like 0{Q~ l ) for fast rotation while 
the x\\ vanishes like 0(Q 3 ). Hence, for isotropic turbulence 
under the influence of a global rotation, the radial eddy heat- 
flux peaks at the poles, and the latitudinal heat-flux goes to the 
poles. 

2.2. Anisotropic turbulence subject to rotation 

The situation is more complicated if the turbulence-field with- 
out rotation is already anisotropic. The anisotropy (unit) vector 
may be G. In this case extra terms appear in which in the 
simplest case run with the second order in G, i.e. 



Xu = ... + X i(GQ)%j +X2(G£i)(GiQj + GjOd. 
These terms provide the new contributions 

Xrr = ■ • ■ + Of! + 2X2) COS 2 9Q 2 



X0r 



- X2 sin 6 cos 6Q 



(8) 



(9) 



2.3. Rotating free turbulence 

There is a close relation of the heat-flux tensor and the one- 
point correlation tensor 



Qij = <«;(*, f)w}(*,0>- 



(11) 



In order to demonstrate this basic outline we start from the 
quasilinear relation 



XU - J 



I 



-io> + xk 2 
Xk 2 Qijdkdu> 
co 2 + x 2 k 4 



dkdcu = 



+ i 



u)Qijdkdu> 
co 2 + x 2 ^ 



(12) 



between the spectral tensor Qij of the turbulence and the heat- 
conductivity tensor Xij ( e -g- Riidiger 1989). The last term in 
dl2l provides an antisymmetric part of the ^ (/ -tensor which is 
still ignored. 

In the remaining term for ^ — > a Dirac (5-function is in- 
volved so that in this approximation 



Xi 



* j Q u {k,Q)dk= l - j e,y(0,T)dr. 



(13) 



If this relation is applied and if the r-integral is approximated 
by T con then 



1 



Xi. 



1 cori" \c i J 



Qi. 



(14) 



It makes sense, therefore, to discuss briefly the one-point 
correlation tensor Qij under the influence of rotation. In the 
Appendix the over-all structure of this tensor is given within 
the r-approximation which is adopted in the following. Two 
properties are particularly stressed here: i) the anisotropy of 
the rotating turbulence at the poles and ii) the pole-equator dif- 
ference of the radial turbulence intensity (u' 2 ). 

At the poles we have {uf) = (u'^). For the characterization 
of the anisotropy it is thus enough to write 



{u' r 2 )-{u' 2 )-{u' 2 ) = {u r 2 )-2{u' 2 ) = 



-^ + fl(0 1 +0 2 +20 3 -0'))(4 O)2 > 



(15) 
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Fig. 1. Contributions of the isotropic (full line) and anisotropic 
(dashed line) parts of turbulence into the normalized pole- 
equator difference (I19> of the radial turbulence intensity. The 
full and dashed lines are the functions <f>\\ and <f>' + (p\ . 



with a after Eq. (IA.4t . They will be considered in this Section 
exclusively. The azimuthal heat-flux is discussed in Section[6] 
With the relations JA.11> it follows for fast rotation that 



Here 



(uf) - (u?) - «) » ~ <U '"' 2> - 



(17) 

is the Coriolis number (or inverse Rossby number) of the tur- 
bulence. Hence, for isotropic turbulence (a = 0) 



« 2 > - <«; 2 > + < M ; 2 >. 



(18) 



At the poles the (fast) rotation of an originally isotropic turbu- 
lence results in a dominance of the vertical turbulence compo- 
nents. The effect is reduced for radial-type turbulences (a < 0) 
and it is enhanced for horizontal-type turbulences (a > 0). 

Canuto, Minotti & Schilling (1994) demonstrated for their 
rather general model how the global rotation influences the for- 
mation of anisotropy between the components of the turbulence 
intensity. If in the equatorial plane without rotation both the tur- 
bulent intensities (u' r 2 ) and (u'^) strongly differ then the rotation 
is smoothing the differences and may even completely suppress 
them (Riidiger, Tschape & Kitchatinov 2002). For faster ro- 
tation there is thus a clear tendency for a return-to-isotropy . 
A similar (even crossover) behavior for anisotropic turbulence 
has also been found for his rotating convection-turbulence by 
Chan (2001). 

For the radial turbulence intensity one finds the relation 



- « 2 > 



pole 



= (*H+a(0fl + *i)) (uf )2 ). 



(19) 



eq 



The contributions of the isotropic and anisotropic parts of tur- 
bulence to the pole-equator difference dl9> are shown in Fig. [2 
as functions of the Coriolis number. For radial-dominated 
free turbulence (a < 0) the expression dl9l is always posi- 
tive. Therefore, {u' r 2 } for rotating turbulence should be always 
greater at the poles rather than at the equator. 



3. Basic equations, the model 

3D numerical simulations of compressible, thermal convec- 
tion under the influence of rotation are made with the finite- 
difference, fractional-step code NIRVANA (Ziegler 1998, 
1999) in a small rectangular box defined on a Cartesian grid. It 
can be considered as a small piece of a spherical shell. The do- 
main is placed tangentially at different colatitudes from 9 = 0° 
(at the north pole) to 180° (at the south pole) with a step of 30°, 
and contains a convectively unstable layer, surrounded by sta- 
ble stratified layers with overshooting convection. The height 
of the convection zone d is chosen as the unit length. 

The box rotates around the polar axis from west to east (the 
angular velocity vector £1 points toward the north pole). The 
geometry of the computational domain is (x, y, z) e [0, 8c/] e 
[0,8c/] e [-2d, 0]. This volume is discretized by 100x100x128 
grid points which are uniformly distributed in each coordinate 
direction. 

The governing equations describing thermal convection in 
a rotating stratified medium are 



(16) dpu 



-!— = -V ■ (puu) -VP + V • n + pg - 2pQ x u, 

at 



dpU 



= -V ■ (pUu) -PV-U + Q vis + V • (pC p ^Vr). 



(20) 
(21) 
(22) 



The notation for physical variables is standard (U thermal en- 
ergy density, n the viscous stress tensor and Qvis the viscous 
heating term. 

The equations fl20l >— fl22i are closed through the ideal gas 
equation P = (K/fi)pT. The initial distribution of the physical 
quantities represents a 3 -layer poly trophic stratification. We as- 
sume all quantities to be periodic in the horizontal directions. 
At the bottom (z = -2) and top (z = 0) of the box impermeable 
conditions are imposed for the vertical velocity, while the hor- 
izontal velocities satisfy stress-free boundary conditions. The 
temperature and density are fixed at the top of the domain and 
a constant heat-flux is injected at the bottom. 

The dimensionless parameters Ra, Pr and Ta are used to 
control the simulations. In our calculations Ra = 3 • 10 5 and 
Pr = 0.1, while Ta e {10\ 10 6 }. 



4. Rotation-induced anisotropic turbulence 

We start with a discussion of the basic anisotropy between 
vertical and horizontal turbulence intensities without (Fig. |2| 
left) and with rotation (Fig. |2] middle & right). Without rota- 
tion except in the top layer the turbulence is vertically dom- 
inated (Fig. |2 left). As it must be, both the horizontal inten- 
sities are equal. The same is true for Ta = 10 6 in the polar 
region (Fig. 12 middle). For the equator, however, the vertical- 
horizontal anisotropy is more and more reduced but a new 
anisotropy developes between both the horizontal components 
(Fig. El right). 

The simulations also provide the anisotropics in the over- 
shoot region. This is in particular important for the lower over- 
shoot region in relation to the tachocline discussion. It has 
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Fig. 2. Ratios of the turbulence intensity for the nonrotating and rotating turbulence fields. Solid: (u^)/(u'^), dashed: (u'q)I{u'}). 
Left: Except the surface layers the turbulence is strongly vertical-dominated. It is even isotropic in the lower overshoot region. 
Middle: Ta = 10 6 , pole. Right: Ta = 10 6 , equator. Under the influence of rotation only the lower half of the box remains vertical- 
dominated. Note that the basic rotation originates a (mild) dominance of the horizontal motions over the radial motions in the 
lower overshoot region. 



been argued that the stability of this zones against convection 
changes the turbulence to the horizontal-type (Spiegel & Zahn 
1992). This is not observable in Fig.[2]for the case without ro- 
tation (left). 

The behavior of the vertical turbulence intensity is even 
more important. For free and anisotropic turbulence we ex- 
pect for vertically-dominated turbulence that the polar values 
exceed the equatorial values (see Appendix). The opposite be- 
havior is shown by the simulations. Fig. [3] reveals that in the 
bulk of the convection box the equator dominates the poles, i.e. 



« 2 > 



> « 2 > 



(23) 



pole 



This unexpected result (see Eq. W9Y ) has important conse- 
quences for the meridional components of the heat-flux. After 
Eq. (1131 the heat-conductivity tensor is proportionate to the 
one-point correlation tensor (multiplied with a correlation 
time) so that it should not be too surprising if also the radial 
heat-flux F r would peak at the equator rather than at the poles. 
The simulations confirm this expectation. 

5. The poloidal heat-flux 

Figure|4](top) shows the depth-profile of the correlation (u' r T') 
in the box for various latitudes. Due to the rotation the values 
differ for poles and equator. The pole-equator difference, how- 
ever, of the radial heat-flux depends on the radius. Except the 
top layer the eddy heat-flux at the equator exceeds the eddy 
heat-flux at the poles. In the top layers where after Fig. 13 the 
turbulence is horizontally-dominated the polar heat-flux domi- 
nates the equatorial one. 

This is a characteristic but unexpected result. It does not 
contradict, however, the findings of Tilgner & Busse (1997, 
their Fig. 8) in which the latitudinal profiles of the radial heat- 
flux differ for differing parameters. For the cases with large 
Pr at the inner boundary the heat-flux dominates indeed at the 
equator; and at the outer boundary the heat-flux dominates in- 
deed at the poles. 
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Fig. 3. The radial turbulence intensity for Ta — 10 6 . Note that 
except the top layer the turbulence at the equator exceeds the 
polar values. Quite a similar result follows from box simula- 
tions without any density stratification (Giesecke, Ziegler & 
Riidiger 2004). 

Also Rieutord et al. (1994, their Fig. 8a) and Kapylii, Korpi 
& Tuominen (2004, their Fig. 7) found similar results. Here 
we are led to the general conclusion that a crossover exists of 
the pole-equator difference of the radial eddy heat-flux almost 
at the same depth where the vertically-dominated turbulence 
changes to a horizontally-dominated turbulence. As we have 
demonstrated with Eq. Jl 3i the behavior of the radial heat-flux 
is a direct reflection of the rotation-influenced radial turbulence 
intensity (u'}). It is shown in Fig.[3]that in the box (except the 
outermost layer) the {u'}) at the equator exceeds the value at 
the poles. 

A similar crossover does not exist for the latitudinal eddy 
heat-flux {u' g T') plotted in Fig. |4] (bottom). This heat-flux van- 
ishes by definition at the poles and the equator. Between these 
extrema the heat flows towards the pole in the convection zone 
(and towards the equator in the lower overshoot region). It is a 
consequence of the Coriolis force which is not involved in the 
description (0. The function x\\ in dlOi o. obviously, is positive- 
definite as was predicted in Section l2~TT bv a much simpler con- 
sideration. 
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Fig. 4. The correlations {u'T') for different colatitudes in the 
box simulations after horizontally and time-averaging vs depth, 
Ta = 10 6 . Top: «7">, bottom: (u' g T'). Solid: pole, dashed: 30°, 
dotted: 60°, triple-dot-dashed: equator. 
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Fig. 5. The coefficient function^- in the box for different Taylor 
numbers. Note the characteristic zero close to the top of the 
convective box. 

6. The azimuthal heat-flux 

Obviously, the last term in is linear in Q in opposition to the 
second one which only appears in the second order in Q. The 
expression linear in Q leads to an azimuthal component of the 
heat-flux of 



F<j, = pCpXtrPr®, 



(24) 



where j3 r > in convection zones. The sign of^ r , therefore, 
determines the sign of Fq. One could argue, however, that for 
axisymmetric constellations the F^ is not important because of 
divF = in that case but i) many stars are not axisymmet- 



ric and ii) the linear-in-£? term F$ should be important for the 
comparison of (quasilinear) theory and (nonlinear) simulation. 
If the rather simple analytical theory of Fq cannot be confirmed 
by numerical simulations then the chance is small to under- 
stand the higher-order terms in l|6} which concern the poloidal 
components of the vector F. 

The azimuthal heat-flux results from the relation 



X<t, r = sin 6Qx. 



(25) 



The quantity^ in Eq. d25i has been computed for a simple tur- 
bulence model subject to a basic rotation but without stratifi- 
cation by Kitchatinov, Pipin & Riidiger (1994). The term was 
positive-definite but it proved to be very small. The main rea- 
son for its smallness is that the integral kernel runs with u> 2 so 
that in the r-approximation the term vanishes. 

This is not true if anisotropic turbulence is considered. 
After 1141 we have 



(26) 



As known the zonal cross correlation Q$ r only exists for rotat- 
ing stratified turbulence and we write 



* r = vj£3Vsin( 



(27) 



with V as the vertical component of the A-effect. Along this 
way we arrive at 

Xtr=^v T Vsm6. (28) 

The zonal heat-flux 1241 should thus be positive (negative) for 
positive (negative) V. A very similar relation has been formu- 
lated by Gough (1978). Obviously, the F^ and the radial A- 
effect have the same sign. In Fig.|6]the radial A-effect quantity 
V is given obtained by Kitchatinov & Riidiger (1993) within 
the quasilinear r-approximation. Generally, the function V is 
negative which result has recently been confirmed numerically 
by Chan (2001) and Kapyla, Korpi & Tuominen (2004). Hence, 
the F,/, should be negative and this is also the result of the sim- 
ulations shown in Fig. [7] Positive (negative) fluctuations u'. are 
thus be correlated with negative (positive) temperature fluctua- 
tions which should be observable. 



7. Rotation-law consequences 

As mentioned in Section^the rotation laws which have been 
computed so far (Kitchatinov & Riidiger 1999; Kiiker & Stix 
2001) were obtained with a radial heat- flux which peaks at the 
poles as the result of the basic rotation. The simulations do not 
confirm this latitudinal profile. On the other hand, the rotation 
laws are computed with a latitudinal heat-flux which goes to the 
poles and this is confirmed by the simulations. We have thus to 
check the consequences of the new situation for the rotation 
law theory. 

As a demonstration of the complex character of the result- 
ing mean-field rotation laws in Fig. [8] (top) the rotation law is 
given as due to the A-effect alone (details given by Kiiker & 
Stix 2001) without both meridional flow and rotation-induced 
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Fig. 6. The A-effect quantity V as a function of latitude and 
Coriolis number (i.e. Q). It is mainly negative between pole 
and equator. The turbulence model is highly density-stratified. 
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Fig. 7. The correlation (u/T') for different colatitudes in the 
box simulations after horizontally and time-averaging vs depth, 
Ta = 10 6 . 

eddy heat-flux. The rotation profile well complies with the ob- 
servations (see Kiiker, Riidiger & Kitchatinov 1993). 

The inclusion of the meridional flow (Fig. [8] bottom) drasti- 
cally changes the situation. The resulting poleward surface flow 
strongly reduces the equator-pole difference and according to 
the Taylor-Proudman theorem the S-isolines become parallel 
to the rotation axis. In this case one finds a slight superrota- 
tion beneath the equator and a rather uniform angular velocity 
beneath the pole - in great contrast to the observations. 

The situation changes with the rotation-induced eddy heat- 
flux included. The poles become warm and a circulation de- 
velopes towards the equator. Hence, the total meridional flow 
becomes rather slow and again we have the situation of Fig. [H] 
(top). This is the solution of the 'Taylor number puzzle' which 
consists in the existence of two opposite directions of the 
meridional circulation (Fig. [9] top). In Fig. [9] (bottom) the nu- 
merical experiment with (artificially) xer = reveals that it 
is the ^0 r -effect which solves the Taylor number puzzle rather 
than the latitude-profile of^-,-,-. 

8. Discussion 

For rotating free anisotropic turbulence Eq. dl9l provides a 
dominance of the vertical turbulence intensity at the poles. 
Our box simulations reveal the opposite behavior. As shown 
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Fig. 8. The rotation law in the solar convection zone without 
the influence of the rotation-modified heat-flux without (top) 
and with (bottom) meridional circulation. 

in Fig. 13 the turbulence at the equator dominates. This is not 
a consequence of the density stratification which was not in- 
volved in the evaluation of the correlation tensor for rotating 
free turbulence as given in the Appendix. Giesecke, Ziegler & 
Riidiger (2004) considered boxes without any density stratifi- 
cation but with the same result. One needs global simulations 
in order to find out whether or not the limited size of the boxes 
is responsible for the unexpected effect. 

An immediate consequence of this result concerns the 
latitude-dependence of the vertical turbulent heat-transport 
which is given in Fig. |4] (top). In the main bulk of the box 
the heat-flux at the equator exceeds the vertical heat-flux at the 
poles 1 . For free turbulence the quasilinear theory on the basis 
of the r-approximation leads to the opposite behavior. There 
is possibly a simple explanation of the differences for the lati- 
tudinal profile of (u' r 2 ) for free turbulence and for box simula- 
tions. For rapid rotation the Taylor-Proudman theorem strongly 
damps the M r -component of the turbulence at the poles due to 
the boundary condition u' r — at the top and bottom of the box. 

This argument hardly holds for the cross-correlations 
(u' g u' r ) and (u'^u' r ) which identically vanish at the poles (see 
Eqs. 1A. 121 and JA.13» . The results of the simulations for the 
eddy heat-flux components Fq and F$ do fully comply with the 

1 In the top layer of the box where after Fig. [5] the polar value of 
(ii' r 2 ) dominates the equatorial one also the polar value of the heat-flux 
dominates its equatorial value 
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Fig. 9. The same as in Fig.[8]but with the full rotation-modified 
heat-flux (top) and xbt - ( see Eq. ( II Oh bottom). Obviously, 
the rotation-induced off-diagonal component xer produces the 
main effects. 



analytical quasilinear expressions. Both correlations are nega- 
tive at the northern hemisphere 2 . The latitudinal heat transport 
is therefore always polewards and the azimuthal heat-flux is 
always westwards. The results of Kapyla, Korpi & Tuominen 
(2004) for the latitudinal heat-transport are here less coherent. 

The correlation (u'+T') is due to the rotational influence 
and it represents the radial angular momentum transport by 
the turbulence. It should be observable at the solar surface. 
Convective-originated temperature fluctuations like those of 
the mesogranulation are expected to be anticorrelated with their 
local angular velocity fluctuations (see Duvall & Gizon 2000) . 

For axisymmetric stars, of course, the longitudinal heat-flux 
is not important because of dF,p /d(p — but this might not to 
be true for close binaries or also for single stars with distinct 
nonaxisymmetric surface structures like the (flip-flop) FK Com 
stars. 

More clear is the meaning of the latitudinal heat-flux com- 
ponent Fg. It is always polewards directed leading to the for- 
mation of warm polar caps - not in agreement with the results 
of Rieutord et al. (1994). We have shown with the simulations 
presented by the Figs. [8] and [9] how important this effect for 
the explanation of the solar/stellar differential rotation is. The 
meridional flow driven by the warm poles breaks the Taylor- 
Proudman theorem which would otherwise produce cylindri- 



cal isolines of the angular velocity which are not found by the 
the inversions of the helioseismology though. As we have also 
shown, the role of the radial eddy heat-flux F r is here of minor 
importance. 
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Appendix A: Rotating anisotropic turbulence 

The usual approach to rotating turbulences is to prescribe the 
turbulence for nonrotating fluid and then to derive the influence 
of rotation on the given original turbulence. As the first step in 
this procedure, we consider an incompressible fluctuating ve- 
locity field, u (0 \ in a nonrotating fluid. The statistical properties 
of the original turbulence are given by its spectral tensor, 

<® = ^K lj( k) 



F<s is also negative at the southern hemisphere 



I6nk 2 
E\(k, to) 

1671k 4 



(KijiGWUG) - K im (G)K jn (G)) (A.l) 



where G is the radial unit vector, k is the wave vector, to is 
frequency, and K/j is the projection tensor, 

,2 
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+ <pz(Q,*)GiGj - 



K u (k) = Sij - kikj/k' 



(A.2) 



Without rotation, the only preferred direction can be the ra- 
dial one. Equation (IA.U allows for the radial anisotropy. 
Characteristic velocities are given by 



{uf )2 ) = ij E(k,to)dkdto, 
o 

(m^ 2 ) - 2{uf )2 ) = i j Ei(k, to) dkdto, 



(A3) 



where u^ 1 is the horizontal velocity, {u^ 2 ) = {uv 12 ) + (uf )2 ). 



Spectral functions must satisfy E\ 



-E. A dimensionless 



anisotropy parameter, a, is introduced 

a = ((«[ i 0)2 )-2(«r))/<«W 2 ), (A.4) 

with a > -1 (Bochner's theorem). Negative a means anisotropy 
of radial type with (m^ 2 ) = (u g 0)2 ) < (uf ) while positive a 
means predominance of horizontal motions. 

The rotational influence is involved via the tensor 



DuiQ) 



6ij + 



(2k-i>) 



1 + 



(2k-af- 



(A.5) 



k 2 (— ioi+v^k 2 ) 2 



(Riidiger 1989), which expresses the rotating turbulence in 
terms of the original one by the linear relation 



Q u (Jfc, to) = D im (Q) D% {£}) QZ (k, to) , 



where the star here means complex conjugate and v t is the ef- 
fective viscosity by microscale turbulence. 

With JA.6> . the parameters of the rotating turbulence can 
be derived. We make, however, further simplifications by in- 
troducing the mixing-length approximation via a special form 
of the spectral functions 



E (k, CO) = 6 ■ {uf )2 )5 (k - Crr) 8 M . 



E\ - aE, 



(A.7) 



with ( 2 on /v t 



Tcorr (Kitchatinov 1991), where £ corr and r con - 
are the mixing length and time respectively. An advantage of 
iA.l\ is that the effects of rotation can be expressed in terms 
of relatively simple parameters like the Coriolis number. The 
tensor iA.5\ now reads 



Du(Q) 



Sij + cos 9Q.* etjpkp/k 
1 + cos 2 9Q.* 2 



(A.8) 



where 9 is the angle between the angular velocity and the wave 
vector. 

The relation \A.6\ can now be integrated over the wave 
number space to find the one-point-correlation tensor. It reads 



. • G) i x 
i 3 (O ) K —^- [QiGj + a p) 



(tp[ (Q* ) 12 + tf> 3 («')) ((12 ■ G) 2 Sij + Q t Qj) /& 



+ al (O ) [Gi€ jmp + Gj£ imp j -7tG p + 



, (Q-G) / 



a 



jmp Qj^imp 



^ £2 m G p 



for rotating anisotropic turbulence. For slow rotation 
a 1, <p' =* 1/2, tf> 2 =* -1/2, I =* 2Q*/5, 



(A.9) 



(A. 10) 



(Q* «; 1, other functions are of second or of higher order in 
Q*), and for fast rotation (Q* » 1) 



3tt 



9tt 
64Q* 



/o 



3tt 

" 640* 
3tt 



16Q* 



- <P2 

16£l* 



3n 
32Q* 



(A.ll) 



(£2* » 1). Full expressions for the functions are too bulky 
to reproduce here. Some of them can be found elsewhere 
(Kitchatinov, Pipin & Riidiger 1994; Kitchatinov 2004). 
For the off-diagonal component, Qg r yields 



Q 9r = - sinfl cos 9 {uf )2 ){<p\\ (£T) 

+a (ft*) 12 + cos 2 9 4>i (Q*)) ). 



(A. 12) 



Qg r is thus always negative. 

The sign of the off-diagonal component of the correla- 
tion tensor is only controlled by the turbulence anisotropy, i.e. 



(A.6) = {uf )2 )a sin0(/ o (Q*) + cos 2 9 h (O*)) . 



(A.13) 



Q^, r is thus negative for radial-type of anisotropy (a < 0). It 
vanishes for isotropic turbulence (a = 0) and is one of the 
known components of the A-effect. 

With the mixing-length approximation iA.7\ . the relation 
between the turbulent thermal conductivity and the velocity 
correlation tensor becomes 



Xij 



rQij 



(k,to) dkdto. 



(A. 14) 



Note the factor 2 by which this relation differs from dl4> which 
only holds in the limit x — * 0. 



Qij = <4 O)2 )-jl0(Q*)<5, 7 + 



QiQ; 



+ aU'(a*)Sij + tpi (Q*) 



(£1 GY 



